Scaling of critical connectivity of mobile ad hoc communication networks 
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In this paper, critical global connectivity of mobile ad hoc communication networks (MAHCN) is 
investigated. We model the two-dimensional plane on which nodes move randomly with a triangular 
lattice. Demanding the best communication of the network, we account the global connectivity r\ 
as a function of occupancy a of sites in the lattice by mobile nodes. Critical phenomena of the 
connectivity for different transmission ranges r are revealed by numerical simulations, and these 
results fit well to the analysis based on the assumption of homogeneous mixing . Scaling behavior 
of the connectivity is found as n ~ f{R l3 a), where R = (r — ro)/ro, ro is the length unit of the 
triangular lattice and /3 is the scaling index in the universal function fix). The model serves as 
a sort of site percolation on dynamic complex networks relative to geometric distance. Moreover, 
near each critical er c (r) corresponding to certain transmission range r, there exists a cut-off degree 
fc c below which the clustering coefficient of such self-organized networks keeps a constant while the 
averaged nearest neighbor degree exhibits a unique linear variation with the degree k, which may 
be useful to the designation of real MAHCN. 

PACS numbers: 89.75.Hc, 89.20.Hh 
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Mobile ad hoc communication network (MAHCN) [l|, Q 
is a new sort of communication circumstance. It con- 
sists of many mobile nodes carrying out collective duty 
while nodes communicate with each other via wireless 
links. Neither central control authority nor intermediate 
services such as base stations for cellular mobile-phones 
exist in the network. Each of its nodes needs to relay 
packets within its limited transmission range for other 
participants in multi-hop edges. MAHCN changes its 
topology with time withoutprior notice since its nodes 
are free to move randomly Q. Therefore, to realize ef- 
fective communication and to do moving jobs, it should 
self-organize into a dynamically stationary network by 
certain local protocols. The study of mobile ad hoc net- 
works has attracted much attention recently due to their 
potential application in battlefield, disaster relief provid- 
ing, outdoor assemblies and other settings with temporal, 
inexpensive usage. Tens of protocols have been proposed 
by designers. However, investigation on property of the 
connectivity as viewed from statistical physics is still in- 
adequate, which motivates the work in the present paper. 

The theory of complex networks can provide pow- 
erful tools to investigate the mobile ad hoc networks. Xie 
et aZ.[|| analyzed the formation of complex networks in- 
volving both geometric distance and topological degree of 
nodes. Sarshar et aZ.@, H[ noticed that while new nodes 
are added to the existing network, other nodes might 
leave the network rapidly and randomly. They presented 
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results about the possible emergence of scale-free struc- 
ture in ad hoc networks. However, they considered only 
static cases instead of analyzing the influence from the 
motion of nodes. Nemeth and Vattay @ studied the giant 
cluster of such networks, and pointed out that the giant 
component size in the percolation could be described by 
a single parameter — the average number of neighbors of 
nodes. On the other hand, models of percolation on net- 
works fiol . [TjT | were often employed to analyze spreading 
processes, especially epidemics with occupancy threshold 
p c showing drastic transitions. 

A communication network may deliver meaningful ser- 
vices only if the network is well connected, or at least has 
a vast subset that is connected. Therefore, one goal of 
studying the ad hoc network is to find out how the net- 
work can maintain its connectivity. Different from previ- 
ous study @ , we demand global connection of all nodes in 
the MAHCN for the best communication, which means 
a stricter case than site percolation j9l [ToL [ill H3. [l3l. [l4| ■ 
For simplicity, we model the two-dimensional plane on 
which nodes move randomly as a triangular lattice with 
N vertices, so that we mimic round transmission ranges 
with discrete hexagons. We define the probability of 
global connection rj as the ensemble average of n/riQ, 
where n is the number of moving nodes globally con- 
nected to the integrated network, and no is the total 
number of them. Critical behavior of order parameter 
rj is found to rely on both transmission range r and the 
occupancy cr(defined as uq/N) of sites (i.e. vertices in 
the triangular lattice). Scaling behavior of the order pa- 
rameter is verified in the form of rj ~ f(R^cr), where 
R is the reduced transmission range, and (3 is the seal- 
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ing index in the universal function f{x). Moreover, at 
critical thresholds of occupancy, individual nodes self- 
organize into complex networks which display particu- 
lar degree distribution p^|. clustering coefficient flr| and 
averaged nearest neighbor degrees 17| of such dynamic 
communication networks. 

To describe the case of self-organized communication 
of individual nodes, it is assumed that they move on a 
two-dimensional plane in a discrete way, and a simple dy- 
namic ad hoc network model is proposed as follows. On 
the two dimensional triangular lattice(see Fig.l) of size 
L, individual nodes are assumed to distribute randomly 
on the sites of it. In our work, the total number of the 
sites are chosen as N = L 2 with L = 200. And, un- 
der periodic boundary condition, we restrain the motion 
of nodes along edges between sites. At the initial time 
step, assign uq nodes to the sites in the triangular lattice 
randomly. Every site of such a lattice can be occupied 
by only one node or nothing. Considering the dynamic 
topology of the ad hoc network, we suppose that, at every 
time step, each node can move randomly in one of the six 
directions to its neighbor site if it is not occupied. Two 
nodes in the ad hoc network can communicate with each 
other if the distance between them is less than the mini- 
mum of their two transmission ranges [H, . To simplify 
our study, we assume that all the mobile nodes have the 
same fixed transmission power j20l[2l|. Therefore they all 
have equal transmission range r for valid communication 
of the whole network. The transmission range r is an im- 
portant parameter in the designation of ad hoc networks, 
since it is vital to keep the network globally connected. 
Proper adoption of range r could minimize energy con- 
sumption since transmission power is proportional to the 
square of it. Note that neither self- loop nor multiple edge 
is allowed in the network: (1) A node should not com- 
municate with itself; and (2) technically there is no sense 
to open another communication channel between any two 
nodes if they are already neighbors. The motion of nodes 
at a certain value of occupancy forms different configu- 
rations of site-occupation on the lattice, which serves as 
the ensemble for the calculation on the global connec- 
tivity^. e. the probability of the global connection) and 
other averaged quantities. 



An intuitive way to ensure the best communication by 
multi-hop linking is to increase the occupancy of sites on 
the lattice by mobile nodes, while this changes the energy 
consumption of whole the network. In the framework of 
traditional percolation problems, continuously increasing 
site occupancy will pass threshold of site percolation. To 
our knowledge, the term "site" can be used dually: one 
means a block[22j; the other means a vertex[23T|. Obvi- 
ously our usage belongs to the later one. Two sites are 
neighbors when two occupied vertices are directly linked 
by a common bond (edge) or two occupied blocks have 
a common border. Therefore, an indirect connection be- 
tween any two sites means that a path exists from one 




FIG. 1: (Colour online) Triangular lattice with transmission 
range r = ro and r = 2ro, respectively, where ro is the length 
of every edge. 



to the other neighbor by neighbor. In the present model, 
however, every node is located at the center of its own 
transmission range r = ztq (z is a positive integer, and 
ro is the length of an edge of any minimal triangle). All 
the other nodes at sites inside this circle connect with 
it directly. Taking the circle r = 2ro (blue in Fig.l) 
as an example, nodes 2, 3 and 4 inside the inscribed 
hexagon (orange in Fig.l) of it are direct neighbors of 
node 1. Therefore, neighbors here are determined by 
the transmission range and in the sense of topological 
connection, just as what occurred in complex network 
models [ToL fll| . which distinguishes them from those in 
traditional two-dimensional site percolation. It is a spe- 
cial case of topological correlation valid within certain 
geographic distance 0, [H, [H, US HI- The triangular 
lattice in our model is a beneficial setting to describe 
moving nodes and for discrete calculation of network pa- 
rameters. Moreover, in percolation problem one always 
focuses on the probability for a site to be included in the 
giant component which just extends from a border to its 
opposite one in a lattice with finite size. However, for 
an ad hoc network bearing search, rescue, tracing or pre- 
cise attack, the task may concentrate on a few, even a 
single moving target. It may demand global connection 
of all nodes, which is quite different from the percolation 
problem which leaves many nodes scattering outside the 
giant component. 

The order parameter 77, i.e. the global connectivity of 
the MAHCN, is calculated with burning algorithm |23|. 
The evolution of it should rely on the occupancy a of 
the sites, and 77 enhances when a increases. The ratio of 
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the enhancement depends on the number n of the nodes 
which have been connected into the largest dynamic net- 
work at that time step, and it also depends on the num- 
ber of the disconnected nodes, based on the assumption 
of homogeneous mixing 30] of randomly moving nodes. 
Therefore, we have 



— oc n(n - n) (1) 
da 

where no is the total number of mobile nodes. Using the 
definition of T) and getting the effect of the transmission 
range included, we arrive at 



= ff(r)»;(l - rf) (2) 

where g(r) is the function of transmission range r. This 
equation can be solved with the uniform initial condition 
i](p -> 0) = t] : 

^ = i T\ V0 \ -o(r)a ( 3 ) 

In fig. 2, simulation results for the global connectivity 
as a function of the occupancy of the sites on the trian- 
gular lattice are illustrated. They are in good agreement 
with the analytical result of eq.(3) under the condition 
of dimensionless function g(r) ~ r/ro- Actually, it is 
naturally expected by dimension analysis on the expo- 
nent in the denominator: g(r) should have no dimension 
since occupancy a is dimensionless. The difference be- 
tween simulation and analytical results at bottom parts 
can be attributed to the deviation from homogenous as- 
sumption by the distribution of nodes at discrete sites on 
the triangular lattice, and size effect. 



Numerical results display the critical behavior of global 
connectivity, i.e. they show its drastic transitions oc- 
curring at critical values a c for different transmission 
ranges. When occupancy a passes <r c , the dynamically 
moving nodes self-organize from a disconnected state 
to a surely globally connected one. For our triangu- 
lar lattice model, a c = 0.37,0.21,0.13,0.09 and 0.065 
for r = 2ro, 3ro, 4ro, 5ro and 6ro, respectively. Obvi- 
ously, various critical values of node occupancy <r c are 
required to ensure global connection for different trans- 
mission ranges in MAHCN, which means that we can also 
inversely choose proper transmission range to minimize 
energy consumption of the network for different density 
of nodes on the lattice. 

It is natural to rescale ry(cr, r) into a universal scaling 
function from direct observation of Fig. 2. we have 
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FIG. 2: (Colour online) Comparison between theoretical anal- 
ysis and numerical simulation on the probability of global 
connection as a function of a for various transmission ranges 
r. Hereafter all simulation results are averaged over 300 re- 
alizations of configurations formed by moving nodes on the 
triangular lattice. 

with reduced transmission range R — (r — r )/r . the 
scaling index (3 of the universal function /(x), respec- 
tively. We draw transition curves in Fig. 3 to show rescal- 
ing process. We can see calculated curves for all r col- 
lapse into the one with r = 2ro, and the index f3 = —0.49 
gives perfect convergence of all the curves. This provides 
the evidence that the transitions at cr c (r) are really crit- 
ical phenomena. The inset of Fig. 3 shows rj versus a for 
different sizes of MAHCN with r = 2r . The critical 
value <7 C is independent of the sizes of lattices, which is 
also valid for different transmission ranges. Near criti- 
cal points <r c (r), nodes with autonomic communication 
self-organize into time-varying complex networks which 
are reminiscent of directed dynamic small- world network 
(DDSWN) model[3l[ but with different scaling variables. 
Indeed, the ratio of the number of nodes receiving mes- 
sage to the total number of nodes should vary in the same 
way as that model provided the nodes move at the same 
speed, have uniform transmission range and relay mes- 
sage without delay, which will be discussed under another 
title. But in the present work we check global connec- 
tivity with burning algorithm, assuming that "combus- 
tion" (similar to message spreading out) is much faster 
than variation of topological structure. Moreover, the 
feature of transmission range-dependence distinguishes 
itself from DDSWN model. It is also noticeable that 
Hu and Chen [33] investigated scaling functions for bond 
random percolation on honeycomb lattices with differ- 
ent aspect ratios. By comparison, the present model is 
pertaining to maximally connected network checked with 
combustion algorithm on a triangular lattice although we 
always pay attention to the hexagonal cell within the cir- 
cular communication range. 
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When the occupancy of nodes just exceeds the value 
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FIG. 3: (Colour online) Scaling behavior of critical global con- 
nectivity of mobile nodes for different transmission ranges. 
77(ln<r) curves collapse into the same one with r = 2ro. Inset: 
Critical connectivity under r — 2ro with different sizes of the 
lattice. 

of a c for certain transmission range r, they are found to 
self-organize into a dynamically stationary network with 
our simulations. We can characterize the connection of 
an ad hoc network with parameters of complex networks. 
The simplest and the most intensively studied parameter 
is degree distribution p(fc)[l5[ because it governs funda- 
mental properties of the system. Degree k of a node, as 
well known, is the total number of its topological edges 
connecting with others. The dispersion of node degree 
is characterized by the distribution function p(k) which 
gives the probability that a randomly selected node has 
exactly k edges. In the present paper, we study param- 
eters of the network when it consists of almost all the 
nodes of MAHCN together. The degree distribution p(k) 
follows Poisson distribution which is different from that 
in ref.(fj, [j], |8|. Fig. 4 shows the ensamble averaged de- 
gree distributions for all simulated transmission ranges r 
when the global connectivity of nodes is above 0.9995. A 
cut-off degree k c = 17 appears in it, which means that 
the probability for any node to have degree k > k c is very 
low. The averaged degree of the whole network can be 
obtained from direct observation of Fig.l, that is, 

< k>= 3r{r + l)cr (5) 

where r = zro as mentioned above. We recognize 
MAHCN from this kind of Poisson-like distribution as 
random networks or small- world networks. By the way, 
the global connectivity rj is also a single- variable function 
of average degree < k > for certain r, which is qualita- 
tively like the behavior of component size S in ref. [9j since 
a is proportional to < k > in such cases. 



Clustering coefficient ( '[jjj] and fc„„[I3; the averaged 
nearest neighbor degree of nodes depicts complex net- 



FIG. 4: (Colour online) The degree distribution p(k) of 
MAHCN for different transmission ranges r. Inset: Cut-off 
degree k c versus varying occupancy a for r — 4ro. 



works as viewed from correlations. For a node i, cluster- 
ing coefficient Ci can be defined as the fraction of pairs of 
node z's neighbors that are also neighbors of each other in 
the topological sense. C(k) of the network is the cluster- 
ing coefficient averaged over nodes with the same value 
of degree k. In our case, C(k) = 0.52,0.55,0.56,0.57, 
and 0.57 for transmission range r = 2ro, 3ro, 4ro, 5ro and 
6ro, respectively, and keep invariant for k < k c . The 
fc-independent behavior of C(k) can be understood from 
the symmetry of hexagonal cells and homogeneous dis- 
tribution. Let us scrutinize the hexagonal cell (orange 
in Fig.l) for r = 2ro as an example. Site 1 in it has 
18 neighbors (assuming full occupation). Therefore, the 
largest possible number of closed topological triangles 
in the sense of complex network should be Cf g which 
serves the denominator of the clustering coefficient of 
it. Linking occupied sites 1 and 2, we search for the 
third one with the distance less than 2tq to both of them 
in the hexagonal cell, which makes 7 triangles. Link- 
ing occupied sites 3 and 4 to the center 1, makes 8 and 
12 triangles, respectively. Therefore, all equivalent sites 
((7 + 8 + 12) x 6/2) in the cell make 81 closed triplets un- 
der the constraint of maximum distance 2ro, and yields 
the clustering coefficient of the site as 81/Cf s = 0.53. 
Under the assumption of homogeneous mixing, this is 
also valid for any occupancy or averaged < k > since 
we only need to multiply both the numerator and the 
denominator of it by local a simultaneously, so that we 
have the same C(k) for small k(k < k c ). The horizontal 
line for k < 17 (see Fig. 5) of C(k) indicates particu- 
lar constant clustering coefficient of MAHCN for degrees 
occurring in high possibilities. Here value k c reflects the 
limit case of connection, and it may be pertinent to the 
structure of the hexagonal cell imbedded in the triangu- 
lar lattice. The high value(above 0.5) of C(k) implies 
that we have dynamic small world networks at critical 
points, and there is large redundancy in communication if 
only strategy of increasing occupancy of sites is adopted. 
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FIG. 5: (Colour online)The clustering coefficient C(k) for dif- 
ferent transmission ranges r on the triangular lattice. 



by the best communication, critical global connectivity 
is found from simulations for various transmission ranges 
by adjusting the node occupancy of sites on the lattice. 
The order parameter scales with transmission ranges, but 
behaves differently from other models. This forms a kind 
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The averaged nearest neighbor degree of node i is simply: 
k n n,i = ^2kj/ki, where j is a neighbor of node i. And 

3 

k n n{k) can be accounted as the function of degree k in 
the following form: k nn (k) = ^ k nn _i/ ^ 1, where V 

fci£V ki£V 

is the subset of nodes with the same degree k. From fig- 
ure 6 we can see that the curve of k nn versus k suggests 
an empirical formula in the linear form for our MAHCN, 
that is 

k nn (k) = b(r) + C(k)k (6) 

where b(r) is a k-independent constant. The positive as- 
sortativity (i.e. increase behavior of k nn (k)) gives the 
particular feature distinguishing it from most other tech- 
nical networks [H, HH, [HI]. The appearance of tails in 
large degrees (k > 17) (see Fig. 5 and Fig. 6) is due to 
very low probabilistic occurrence in the simulation on 
300 realizations. 



Interestingly, the model can be thought as a special 
example for the process of coevolutionary competitive 
exclusion [3^]: taking sites as the state variable of mov- 
ing nodes, each node tends to compete for links to others, 
which is demanded by the global connection. Only nodes 
within other ones' transmission ranges and with the de- 
grees of both ones less than k c can successfully gain links, 
which constitutes threshold conditions of the co-evolution 
network. 

In summary, we present a model for mobile ad hoc 
communication networks by considering uniform trans- 
mission range of nodes and assigning moving nodes ran- 
domly on the plane of the triangular lattice. Demanded 



FIG. 6: (Colour online)The averaged nearest neighbor de- 
gree of nodes k nn (k) versus degree k for different transmission 
ranges r on the triangular lattice. 



of percolation in dynamic complex networks pertaining 
to geographic distance. Moreover, cut-off degree k Ci in- 
variant clustering coefficient and linear assortativity as 
functions of degree are found for self-organized complex 
networks near critical global connectivity. Our model 
suggests that transmission range of nodes and average 
occupancy on the plane should adapt to each other to bal- 
ance minimization of energy consumption with the global 
connectivity. In fact, nodes relaying message consume 
energy continuously. Transmission range of each node 
usually reduces with time relating to its job load. There- 
fore, the present model is applicable only for the rout- 
ing strategy of broadcasting [37|- A more practical model 
should include random distributed nodes with changing 
transmission ranges. Meanwhile, nodes are not necessary 
to move along edges of the triangular lattice. Therefore, 
a random graph model without any lattice is necessary 
for better investigation on mobile ad hoc communication 
networks, which leaves our further work in the future. 
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